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ate  distribution  (MIFR)  if  and  only  if  E[h(x,T)]  is  log  concave 
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§6.  Relation  to  the  MI  FRA  class. 

In  Block  and  Savlts  (1980),  we  defined  a  multivariate  increasing  failure 
rate  average  (MI FRA)  claas  as  follows.  A  nonnegative  random  vector  T  Is 
said  to  be  MIFRA  if  E[h(T)]  <_  E1^<*[h°(T/o)]  for  all  0<a<land  all  non- 
negative  nondecreasing  functions  h;  equivalently,  P(T  e  A)  c^P^^CC e  aA) 
for  all  0  <  a <  1  and  all  upper  sets  A.  It  is  well  known  in  the  univariate 
case  that  if  T  is  IFR,  then  T  is  IFRA  (see  Barlow  and  Proschan  (1975)). 

It  seems  natural  to  expect  the  same  conclusion  in  the  multivariate  setting. 
Although  we  believe  this  to  be  the  case,  we  only  have  some  partial  results. 

Let  T  be  MIFR  and  y  its  Induced  measure.  Let  A  be  an  upper  convex 
set.  Then  for  a > 0,  h(a,t)  -  IA(t/a)  is  log  concave  in  (u,t)  and  non¬ 
decreasing  in  t  for  each  fixed  a.  We  may  assume  without  loss  of  generality 
that  P(T  >  0)  ■  1;  otherwise  POT^-0)  >  0  for  some  i.  But  ^  is  IFR  and  so 
P(Ti  -  0)  »  1.  We  are  thus  reduced  to  the  case  of  one  less  dimension. 

Since  T  is  MIFR,  it  follows  that  G(a)  -  E(h(a,T)]  is  log  concave  in  o>0. 
Also,  since  P(T>0)  -  1,  we  get  that  G(Q+)  -  1.  Consequently,  we  conclude 
that  -  log  G(a)  is  star-shaped  and  so  y(aA)  _>  iia(A). 

This  result  coupled  with  the  fact  that  every  upper  set  can  be  approx¬ 
imated  by  a  finite  union  of  upper  convex  sets  leads  us  to  the  following 
conjecture. 

(6.1)  Conjecture.  If  T  is  MIFR,  then  T  is  MIFRA. 


83  02  014  209 
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Nov  l«t  e  +  0. 

(5.7)  Corollary.  A  nonnegative  random  vector  T  la  MIFR  if  and  only  If 
E(h(x,T)]  la  log  concave  in  x  for  all  functlona  h  which  are  log  concave  in 
(x,t)  and  are  nondecreasing  in  t  for  each  fixed  x. 
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fixed  X.  Let  h(X,y)  •  exp{-4(X,y) }.  Then  h  is  log  concave  in  (X,y), 
and  continuous  and  nondecreasing  in  y  for  each  fixed  X.  Consequently, 
hn(X,y)  -  hn(X,y)  has  the  saae  properties  for  each  n.  Using  the  fact  that 
T  is  MIFR  and  that  hn(X,y)  -*•  lD<X,y)  ai  n  +  »,  ve  obtain  the  inequality 

(2.4) . 

(5.5)  Lemma.  Let  y  be  a  probability  measure  and  A  an  upper  convex  set. 

Then  given  e  >  0  there  exists  a  closed  upper  convex  set  C  c  a  such  that 
y(C)  >_y(A)  -  e. 

Proof .  Given  e  >  0,  choose  a  compact  set  Re  A  such  that  y(K)  ^y(A)  -  e. 

If  G  is  the  convex  hull  of  K,  then  G  is  a  compact  convex  set  satisfying 
Kc  Gc  A.  As  in  Block  and  Savits  (1980),  define  C*{y  +  t:  y  e  G,  t^O). 

It  is  not  hard  to  show  that  C  has  the  desired  properties;  i.e.,  C  is  a 
dosed  upper  convex  set  satisfying  KcGcCca.  Consequently  y(C)  ^y(A)  -  e. 

(5.6)  Theorem .  If  y  is  a  probability  measure  satisfying  the  inequality 

u[XA+  (l-X)B]  >_  uX(A)yX  X(B)  for  all  0<X<1  and  all  closed  upper  convex 

sets  A,B,  then  the  inequality  remains  valid  for  all  upper  convex  sets 

/ 

A,B. 

/ 

Proof.  It  suffices  to  consider  upper  convex  sets  A  and  B  having  positive 
y  measure.  Given  e  >  0,  sufficiently  small,  we  choose  closed  upper  convex 
sets  C  c  a,  D c  B  as  in  Lemma  5.5  such  that  y(C)  ^y(A)  -  e,  y(D)  ^y(B)-e. 
Hence 


y[XA+  (l-X)B]  1  y[XC+ (l-X)D] 


>  yX(C)y1_X(D) >  [y(A)-e]X[y(B)-c]1_X. 
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h(Xx  +  (l-X)x’ ,  Xw+  (l-X)w'  ]  ■  p{D[Xx  +  (l-X)x' ,  Xw+(1-X)w']} 

>,  y{XD(x,w)  +  (l-X)D(x’ ,w')  } 

>.  yX[D(x,w)]y1_X[D(x’ ,w')l  -  hX(x,w)h1“X(x' ,w') ; 

i.e.,  h  Is  log  concave  in  (x,w).  Thus  by  Theorem  2.4,  G  Is  log  concave  In 
x. 

We  now  consider  the  converse.  Let  A  and  B  be  upper  convex  sets.  For 
0<^X_<1,  let  ■  XA+(1-X)B.  We  define  D  ■  {(X,y):  O^X^l  and  yeC^}. 

Then  0  Is  convex  and  so  h(X,y)  -  1^(1, y)  Is  log  concave.  It  Is  not  hard 
to  show  that  h  Is  nondecreasing  In  y  for  each  fixed  X.  By  hypothesis, 
then,  g(X)  ■  |  h(X,y)y(dy)  Is  log  concave  In  X.  Hence 

y[XA  +(1-X)B]  -  g(X)  ^.g^ (l)g1_X(0)  ■  yX(A)y1_X(B) . 

(5.3)  Theorem .  Let  T  be  a  nonnegative  random  vector  and  y(dy) ■  P(T  e  dy) 
be  Its  Induced  measure.  Then  T  Is  MIFR  If  and  only  If 

(5.4)  y[XA-Kl-X)B]  >yX(A)y1’X(B) 
for  all  0  <  X  < 1,  all  closed  upper  convex  sets  A,B. 

Proof.  The  same  proof  as  In  Theorem  5.1  shows  that  If  (5.4)  holds  then  T 
is  MIFR;  If  F(x,y)  is  also  continuous  in  y  for  each  fixed  x,  then  the 
sets  {y:  F(x,y)  z)  are  also  closed. 

Now  let  A,B  be  closed  upper  convex  sets.  As  before,  let  ■  XA+ (l-X)B 
for  O^X^l  and  set  D«{(X,y):  O^X^l,  yeC^}.  Now  D  is  a  closed  convex 
set.  We  set  $(X,y)  •  p(X,y) ,D)  where  p  Is  the  metric  p(u,v)«max  l^-vj.  Note 
that  $  Is  continuous  and  convex.  It  is  also  a  nonincreasing  function  in y for  each 


«*»  1 
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15.  An  alternative  condition  . 

In  this  section  (which  is  somewhat  technical)  we  derive  an  alternative 
condition  that  T  be  MIFR;  it  is  directly  expressable  in  terms  of  the  induced 
measure  y(dy)  -  P(Tedy).  This  result,  coupled  with  some  approximation  ideas, 
allows  us  to  remove  the  continuity  assumption  on  h  in  Definition  4.1. 

He  first  recall  that  a  set  AcJln  is  said  to  be  an  upper  set  if  when- 
ever  xeA  and  y ^x,  then  ye  A. 

(5.1)  Theorem.  Let  y  be  a  finite  measure.  In  order  that  |F(x,y)y(dy)  be 

log  concave  in  x  for  all  log  concave  functions  F(x,y)  which  are  nondecreasing 
in  y  for  each  fixed  x  it  is  necessary  and  sufficient  that 

(5.2)  y[XA+(l-X)BJ>yX(A)y1"X(B) 


for  all  0<X<1  and  all  upper  convex  sets  A,B. 

Proof .  We  shall  closely  follow  the  argument  given  in  Brascamp  and  Lieb  (1975) 
Suppose  now  that  y  satisfies  (5.2)  and  let  F(x,y)  satisfy  the  conditions  of 
the  theorem.  For  each  real  z^O,  let  C(x,z)  «  {y:  F(x,y)  >^z}  and  set 
g(x,z)  -  y [C(x,z) ] .  Note  that  C(x,z)  is  upper  and  convex;  furthermore, 

G(x)  -  |  F(x,y)y(dy)  -  j  g(x,z)dz.  If  we  make  the  change  of  variables 


w 

z-e  , 


then  G(x) 


x)  »  |  h(x,w)ewdw  where  h(x,v)  *  g(x,ew)  -  ylC(x,ew)].  We  also  define 


D(x,w)  ■  C(x,e  ).  It  is  easy  to  show  that  D[Xx+ (l-X)x' ,  Xw+ (l-X)w’]  o 
XD(x,v) + (l-X)D(x’ ,w’)  for  all  0<X<1,  x,x',w,w'.  Consequently,  assuming  (5.2) 
we  obtain  the  inequality 


ponentlally  distributed  random  variables  ,  S{ 

{l,...,n}  such  that  T^  -  minj€j  .  But  (S1 , 
4».(s)  -  min  a  (l<i<m)  satisfy  the  hypothesis 

J£Ji  J 

(4.6)  Example  2.  If  T  has  a  density  which  is  log 


and  subsets  J, , . . . , J 
t  x  m 

. . . ,S  )  is  MIFR  and 
n 

of  Theorem  4.3(v). 
concave,  then  T  is  MIFR 


Apply  Theorem  2.4. 


■"  - 


1 


Si  "  1  i».  with  T  MIFR,  then  E[h(x,S)]  -  E[g(x,T)]  is  log  concave 


(ii)  Let  h(x,s,t)  be  log  concave  In  (x,s,t)  and  continuous  and  nondecreasing 
in  (s,t)  for  each  fixed  x.  We  first  assume  that  h  is  bounded.  Since  T  is 
MIFR  it  follows  that  g(x,s)  ■  E[h(x,s,T)]  is  log  concave  in  (x,s).  Also, 

g  is  continuous  and  nondecreasing  in  s  for  each  fixed  x.  Thus,  since  S  is 
MIFR,  E[g(x,S) ]  is  log  concave  in  x.  But,  by  Fubini,  E[g(x,S)]  *  E[h(x,s,T)]. 

If  h  is  not  bounded,  consider  instead  hAn  and  pass  to  the  limit  (see  Theorem 
2.3(D). 

(i)  If  Jc{l,...,n},  let  ij/j(t)«tj  for  j  e  J.  According  to  Theorem  4.3(v), 

(ipj  (T)  *  T ^  ;  j  e  J}  is  MIFR. 

(iii)  Since  S  and  T  are  independent  MIFR,  (S,T)  is  MIFR  by  Theorem  (4.3(ii). 

Now  take  i|/j(s,t)  *  sj  +  Cj  and  use  Theorem  4.3(v). 

(iv)  Use  Theorem  4.3(v)  with  <Jjj(t)  «  a^t^. 

Vv)  Let  h(x,t)  be  a  bounded  log  concave  function  in  (x,t)  which  is  nonde¬ 
creasing  and  continuous  in  t  for  each  fixed  x.  Since  T^  T  in  distri¬ 
bution,  Efh(x,T  )]  -*•  E{h(x,T) }  as  n  +  «.  But  for  eachn,  E[h(x,T  )]  is 
.  _n  ,  -  -  ~n 

log  concave  in  x.  Consequently,  E[h(*,T>]  is  log  concave  in  x.  If  h  is  not 
bounded,  use  the  argument  as  in  the  proof  of  part  (ii). 

(4.4)  Corollary.  If  T^. . . ,T^  are  independent  IFR  random  variables  and 

4»(t, ,...»t  )  is  continuous,  nonnegative,  nondecreasing  and  concave,  then 
l  n 

4»(T. ,  • . .  »T  )  is  IFR. 

X  n 

(4.5)  Example  1.  The  Marshall  and  Olkin  MVE  (1967)  distribution  is  MIFR. 

This  follows  since  if  T  -  (T, «... ,T  )  is  MVE  then  there  exist  independent  ex- 

_  x  tn 


j 
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§4.  A  multivariate  IFR  class. 

(4.1)  Definition.  Let  T  be  a  nonnegative  random  vector.  Ue  say  that  T  has 
a  multivariate  Increasing  failure  rate  (MXFR)  distribution  If  and  only  if 
E[h(x,T)]  Is  log  concave  In  x  for  all  functions  h(x,t)  which  are  log  concave 
in  (x,t)  and  nondecreasing  and  continuous  In  t ^0  for  each  fixed  x  ^  0. 

(4.2)  Remarks . 

(I)  Again  note  that  according  to  Theorem  2.3(111)  we  need  only  consider 
functions  h(x,t)  with  x  a  single  variable  Instead  of  a  vector. 

(II)  In  Section  5  we  shall  show  that  the  continuity  assumption  is  un¬ 
necessary. 

The  class  of  MIFR  distributions  has  many  desirable  closure  properties. 

(4.3)  Theorem. 

(I)  If  T  Is  MIFR,  then  so  are  all  marginals. 

(II)  If  S  and  T  are  independent  MIFR,  then  (S,T)  is  MIFR. 

(ill)  If  S  and  T  are  independent  MIFR  of  the  same  dimension,  then  S+T  is  MIFR. 

(iv)  If  (T. , . . . ,T  )  Is  MIFR  and  a  >  0  (1  -  1,. . . ,n) ,  then  (a,T, . a  T  )  is 

in  I—  linn 

MIFR. 

(v)  If  T  is  MIFR  and  •  »^m  are  continuous,  nonnegative,  nondecreasing 

and  concave  functions,  then  (<J»,  (T)  , . . .  ,<p  (T))  is  MIFR. 

l  «  m 

(vi)  If  T. ,T.,...  are  MIFR  and  T  converges  to  T  is  distribution,  then  T 

.1  ~t-  .n 

is  MIFR. 

Proof . 

(v)  Let  h(x,s)  be  log  concave  in  (x,s)  and  continuous,  and  nondecreasing  in 

s  for  each  fixed  x.  If  are  as  in  (v)  ,  then  g(x,t)  -  h(x,i|>,  (t) , . . .  ,<|>  (t)) 

«  i  in  „  „  „.  i  „  m  «. 

has  the  same  properties  as  h  (see  Theorem  2.3(11)).  Consequently,  if 


Proof .  Let  R(t)  »-  log  F(t)  be  the  hazard  function  of  T.  We  set  a» 

suptt^O:  F(t)  ■  1}  and  b  "  inf{t^0:  F(t)  ■  0}  (inf  0  ■  +  •).  If  a-b, 

we  simply  take  0  =  a  and  we  are  done;  otherwise  0£a<b£-f<».  Since  R  is 

convex  and  finite  on  (-®,b)  it  is  continuous  there;  furthermore,  it  easily 

follows  that  R  is  strictly  increasing  on  [a,b).  Let  A*lim  R(t)  <_+  <*>.  If 

t-*- 

<p  denotes  the  restriction  of  R  to  [a,b),  then  its  inverse  $  is  continuous, 
strictly  increasing  and  concave  on  [0,A).  The  function  0  is  defined  by 
0(8)  *  inftt^O:  R(t)>s}.  Clearly  0(s)  -  <p  ^(s)  for  0  <_s  <  A  and  0(s)  *  b 
for  s  ^  A.  It  is  not  hard  to  show  now  that  0  has  the  desired  properties. 

(3.4)  Theorem.  T  is  IFR  if  and  only  if  E[h(x,T)]  is  log  concave  in  x  for  all 
functions  h(x,t)  which  are  log  concave  in  (x,t)  and  are  nondecreasing 
in  t  for  each  fixed  x  >_  0. 

Proof.  Suppose  T  is  IFR  and  let  h  be  as  in  the  statement  of  the  theorem. 

Then,  according  to  Lemma  3.3,  there  exists  a  continuous  nonnegative  nonde¬ 
creasing  concave  function  0  such  that  0(S)  and  T  have  the  same  distribution, 

where  S  has  the  standard  exponential  distribution.  Hence 

00 

E[h(x,T)]  -  E[h(x,0(S) ) ]  -  j  h(x,0(s))e“Sds 

0 

is  log  concave  in  x.  This  follows  from  Theorems  2.3  (i),(ii)  and  2.4. 

To  prove  the  converse,  let  h(x,t)  «  l(  .(t).  Then  h  is  log  concave  in 
(x,t)  and  nondecreasing  in  t  for  each  fixed  x.  By  assumption,  then, 

F(x)  -  E[h(x,T)]  is  log  concave  in  x.  According  to  Theorem  3.1,  T  is  IFR. 

(3.5)  Remark .  According  to  Theorem  2.3(111),  we  need  only  consider  functions 
h(x,t)  with  x  a  single  variable  instead  of  a  vector. 


4. 
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S3.  A  new  univariate  IFR  characterization. 

The  univariate  class  of  increasing  failure  rate  (IFR)  distributions  has 
played  an  important  role  in  the  mathematical  theory  of  reliability  (see 
Barlow  and  Proschan  (1975)).  We  shall  first  brief ly  review  the  definition. 

Let  T  be  a  nonnegative  random  variable  with  survival  probability 
F(t)  *  P(T  >  t) .  Set  b  *  inf{t^0:  F(t)  -  0}  (inf  0®+  ®).  We  say  that  T 
has  an  IFR  distribution  if  F(s+t)/F(t)  is  nonincreasing  in  t €  [0,b)  for 
each  s^O.  It  is  well  known  that  the  following  conditions  are  equivalent 
(cf.  Barlow  and  Proschan  (1975)). 

(3.1)  Theorem .  The  following  conditions  are  equivalent: 

(i)  T  is  IFR. 

(ii)  F  is  a  Polya  frequency  function  of  order  two  (PF2) •  i.e.,  F^O  and 

F(x.-y.)  F(x.-y-) 

L  0 

F(x2-yi)  F(x2-y2) 

for  all  -  ®  <  x2  <  ®,  -®<  y^  <  y2  <  ®. 

(iii)  F  is  log  concave. 

(3.2)  Remark.  If  F  has  a  density  f,  then  T  is  IFR  if  and  only  if  the  hazard 
rate  r(t)  ■  f(t)/F(t)  is  nondecreasing  on  [0,b). 

Before  we  state  our  new  characterization,  we  need  the  following  simple  result. 

(3.3)  Lemma.  If  T  is  IFR,  then  there  exists  a  continuous  nonnegative  non¬ 
decreasing  concave  function  \p  on  [0,®)  such  that  <|>(S)  has  the  same 
distribution  as  T,  where  S  is  distributed  as  a  standard  exponential. 
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The  following  facts  about  log  concave  functions  are  easily  verifiable. 

(2.3)  Theorem . 

(i)  If  f  and  g  are  log  concave,  so  are  f  *g,  fAg,  cf  and  f*  for  all 
a  >  0,  c ^0. 

(ii)  If  f  is  log  concave  and  nondecreasing  and  <|t  is  concave  (and  nonnegative) 
then  the  composition  f«tj)  is  log  concave. 

(ill)  f  is  log  concave  if  and  only  if  for  every  x,y  >_  0,  the  function 
g(t)  *  f[tx  +  (l-t)y]  is  log  concave  on  0  <_  t  <_  1. 

(iv)  If  f  is  log  concave,  then  for  every  z,  the  sets  (f  >  z)  and  {f^z} 
are  convex. 

(v)  If  f  is  log  concave,  then  f  is  continuous  on  the  interior  of  the  set 
(f  >  0}. 

We  close  this  section  with  a  very  important  result  about  log  concave  func¬ 
tions.  This  result  is  sometimes  known  as  the  Prdkopa  Theorem  (1971). 

An  independent  and  simpler  proof  is  given  in  Brascamp  and  Lieb  (1975)  . 

» 

(2.4)  Theorem.  Let  F  be  log  concave  on  Rm  x  Rn.  Then  G(x)  *  F(x,y)dy 

is  log  concave  on  Rm.  (Here  dy  is  Lebesgue  measure  on  Rn) . 

The  above  plays  a  crucial  role  in  our  development. 
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§2.  Los  concave  functions. 


Let  A  be  a  convex  aet  in  It  and  £  a  nonnegative  function  defined  on 


A.  We  say  that  f  is  log  concave  on  A  if 


(2.1) 


f[Xx+(l-X)y]  >  fX(x)f1"X(y) 


-Q(x) 

for  all  0<X<1,  all  x,y  e  A.  Sometimes  we  write  f(x)»e  '  where  Q(x)  is 
convex,  but  with  the  understanding  that  Q  may  assume  the  value  +  ®. 

(2.2)  Examples 

(i)  If  Q(x)  is  twice  continuously  differentiable  on  an  open  convex  set 

A  and  has  the  property  that  at  each  point  x  e  A,  the  matrix 

32Q(x)  '  -Q(x) 

[ - — ]  is  nonnegative  definite,  then  f(x)  ■  e  'is  log  concave 

aV*j 

on  A.  Hence  all  Gaussian  densities  are  log  concave. 

(ii)  If  A  is  any  convex  set  in  Rn  ,  then  the  indicator  function  I  (x) 

A  •> 

is  log  concave  on  Hn. 


(ill)  If  f  is  a  nonnegative  concave  function,  then  f  is  log  concave. 


It  is  convenient  to  make  the  following  simple  observation.  Let  A  and 
B  be  convex  sets  in  Kn  and  suppose  that  f  is  log  concave  on  A.  Clearly 
then,  f  is  log  concave  on  B  if  Be  A.  On  the  other  hand  suppose  AcB.  Then 
the  function  f,  which  is  defined  to  be  equal  to  f  on  A  and  zero  on  B\A,  is 
log  concave  on  B.  Thus  without  loss  of  generality  we  may  assume,  that  all 


log  concave  functions  are  defined  on  the  same  convex  set.  In  the  context 

of  this  paper,  it  is  natural  to  work  with  the  convex  set  R?  .  Hence,  unless 
.  .  ,*•*?'  * 

otherjrisd  bpecified,  the  term  log  concave  means  log  concave  on  . 
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SI.  Introduction. 


2)  Various  univariate  classes  of  life  distributions  have  been  introduced 
in  the  mathematical  theory  of  reliability. (see  Barlow  and  Proschan  (1975) 

i - - J 

for  a  detailed  discussion);  Recently  there  has  been  much  interest  in  ob¬ 
taining  multivariate  versions  of  these  classes.  Although  there  have  been 

frttn'i  lartS 

many  different  approaches, '(see,-  e.g.,  Cftfcfgvtew  paper  by  Dlock  and  Savits  < 
•(1981),— we  mentien  only  two  because  of  their  nice  closure  properties: 
the  multivariate  IFRA  class  of  Block  and  Savits  (1980)  and  the  multivariate 
NBU  class  of  Marshall  and  Shaked  ( 1982) .  Both  of  these  classes  are  closed 

under  deletion,  conjunction,  convolution  and  weak  limits. 

^  ''O' 

— In  dlls  paper  we-introduee -  a  multivariate  IFR  class  that  has  similar 
closure  properties.  Its  spirit  also  closely  parallels  that “‘f ound  in 
the  two  previously  mentioned  papers.  Other  definitions  of  a  multivariate 
IFR  class  have  been  proposed  by  Marshall  (1975) ,  but  none  possess  all 
desirable  closure  properties. 

4-  _ 
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?  In-  Section  2  We  present  some  preliminary  facts  about  log  concave 
functions.  Section  3  contains  a  new  characterization  of  the  (univariate) 
IFR  class.  The  multivariate  generalization  and  properties  thereof  are 

given  in  Section  4.  A  useful  alternative  condition  is-  delineated 

>z  th* 

in  Section  5.  Finally,  we- 'compare • this  class. with  the  multivariate 

r  f\ 

IFRA  class  of  Block  and  Savits  (1980)  in  Section  6.  ^  “ 


All  functions  and  sets  in  this  paper  are  assumed  to  be  Borel  measureable 
with  respect  to  Rn.  In  most  cases,  however,  we  do  not  specify  the  dimension 


? 
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